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Abstract: We investigate the evolution of the mutual information between two spatial subsys¬ 
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holographic process, given by the 2+1 dimensional Vaidya-BTZ spacetime in global coordinates, 
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mutual information observed in this case. We also consider the analogous process in the adiabatic 
limit and compare these two cases from a geometrical point of view. 
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1 Introduction 

There has always been a wide interest in non equilibrium phenomena in quantum field theory, since 
they can be used to describe processes in various areas of physics, including condensed matter 
and high energy physics. In particular, it is very important to gain a deeper understanding of 
thermalization, meaning the way in which physical quantities attain their equilibrium values after 
some kind of perturbation has acted upon an initial equilibrium state. Of course, QFT implies that 
a pure state will remain pure under time evolution, so it is interesting to understand precisely how 
the final pure state will resemble a thermal one. 

A simple model of a thermalization process is the so-called “quantum quench”, which involves 
a change in the parameters of the Lagrangian of the system, either instantaneously or over a finite 
period of time. For example, it can describe a uniform or point-like energy injection, depending on 
the kind of sources that are turned on. So, the vacuum state of the initial Lagrangian will be a 
high energy state of the final Lagrangian, and it will evolve with the final Hamiltonian. 

It would be desirable to study the above problems in the strong coupling regime of QFTs, where 
perturbation analysis becomes inapplicable. Although some progress had been made using QFT 
techniques, this problem became more tractable with the appearance of the proposal that certain 
strongly coupled d-dimensional quantum field theories with a large number of degrees of freedom 
are dual to classical gravity theories in d -|- 1 dimensions (this idea originated in [l]-[3]; see [4] for 
a recent review). Specifically, this so-called “AdS/CFT correspondence” suggests that the vacuum 
state of a conformal field theory in d dimensions is dual to pure AdS^+i, while a thermal state is 
dual to an (asymptotically AdS^+i) black hole geometry. So, it makes sense that the process of 
gravitational collapse ending in a black hole in an asymptotically AdSd+i spacetime geometry is 
dual to thermalization in a CFT. 
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Apart from various correlation functions or physical observables, a particularly interesting ob¬ 
ject to study is the entanglement entropy. Even though it is not a physical observable and thus not 
directly measurable, it encodes the way quantum information is distributed in the system, as well 
as the quantum correlations between the subsystems. It is defined as the von Neumann entropy 
Si = —tr(/0xln/9i) of the reduced density matrix px corresponding to a subsystem T. Usually we 
take X to be a spatial subregion of our theory and then Si measures the correlations of X with the 
rest of the system. Si is typically divergent but, after regularization (usually by introducing a UV 
cut-off), we can isolate the finite part (this was investigated in [5]; see [6] for a review). One can 
also define the mutual information of two regions Xi and X 2 by I (Xi, X 2 ) = -|- Sx^ — Sx^ijx^. This 

object is interesting in its own right because it quantifies the correlations shared between Xi and 
X 2 but not the rest of the system, it is an upper bound for bounded correlators between Xi and X 2 , 
and it is finite as long as the two subsystems are not adjacent. 

In general, the entanglement entropy is very hard to compute with QFT techniques, and it 
has only been done for special theories (e.g. free scalar field theories, CFTs), in special states (e.g. 
vacuum or thermal) or for specific regions with some degree of symmetry (see for example [7], [8] 
among others). However, it was realized that holography greatly simplifies this problem. 

For generic theories with non-static Einstein-Hilbert gravity duals, the holographic entangle¬ 
ment entropy (HEE) Sx of a boundary region X (not necessarily on a constant t-slice) is computed 
using the “HRT” proposal of [9]. According to this proposal, to order C>(G)^^), the HEE is given by 
the area (in Planck units) of an extremal surface A anchored on the boundary dT of X (i.e. with 
dA = dl) and homologous to X. Specifically, the homology constraint requires the existence of an 
everywhere spacelike manifold with boundary A U X^. In case there are many surfaces satisfying 
these conditions, we should choose the minimal one, so we can write Sx = min^ easy 

to see that if the spacetime is static, the HRT proposal coincides with the “RT” proposal of [12] and 
[13], replacing the extremal surface by a minimal one on a Euclidean constant-time bulk slice. Note 
that the choice of quantum state in which we compute the entanglement entropy is implicit in the 
choice of a bulk gravity solution, and that the regularization procedure analogous to introducing a 
UV cut-off in field theory is the introduction of a large radial cut-off Too in the bulk spacetime. 

In this paper, we are going to investigate the evolution of the holographic mutual information 
(HMI) between two regions in a compact 1+1 dimensional GET after a kind of quantum quench 
describing an instantaneous, uniform energy injection into the ground state of the system. We are 
going to use the simplest possible dual gravitational analogue: the collapse of a spherically sym¬ 
metric, infinitely thin null shell in pure 2+1 dimensional AdS spacetime, resulting in the creation 
of a BTZ black hole. We will ignore the details of the matter which gives rise to the so-called 
global Vaidya-BTZ metric^ describing the above spacetime, focusing only on the purely gravita¬ 
tional part, which is expected to give general and universal results. We will only consider this low 
dimensional model because it is easier to apply the HRT prescription in this case: the extremal 
surfaces are spacelike geodesics. However, it is able to capture important information about the 
higher dimensional cases too. 

Some aspects of our general setup can be seen in figure 1 (left). A spatial slice of pure AdSa 
is presented, and two boundary intervals Xi and X 2 have been singled out. Obviously, the HEE 
Sxi of the interval Xi will be computed by the length of the spacelike geodesic {AB) and Sx^ by 
the length of {CD). However, for the union Xi UX 2 there are more candidates: both the set (AB) 
and {CD) (blue geodesics - “connected” configuration) and the set {BC) and {DA) (red geodesics 
- “disconnected” configuration) satisfy the boundary and homology conditions of HRT, so Sxjui 2 
will be computed by the set with the smallest sum of lengths. Thus, we can write the HMI in the 

^See [10] for a discussion of some subtleties (which imply that this condition requires some refinement). 

^Another model for holographic quenches was proposed in [14], but it will not be considered here. 
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Figure 1. Spatial projection of spacelike geodesics, with the boundary brought to finite distance using 
the coordinate p = tanr. Left: Spacelike geodesics in pure AdSa spacetime, anchored on the endpoints A 
and B of the boundary interval Xi and C and D of I 2 in various combinations. Right: Spacelike geodesics 
in BTZ spacetime with r+ = 1/4. The inner black circle marks the event horizon of the black hole. The 
two blue geodesics wind around the black hole 0 and 1 times, while the two red ones wind around the black 
hole 0 and 2 times. The innermost blue geodesic is homologous to the smaller interval I and the outermost 
red geodesic is homologous to the complementary interval I. 


form 

I(Ii,l 2 ) = {AB) + {CD)-min[{AB) + {CD),iBC) + iDA)] 

= max[{AB) + (CD) - {BC) - iDA),0]. 

We can immediately see that eq. (1.1) implies > 0. It is important to note that the dashed 

geodesics (AC) and (BD) do not satisfy the homology condition (there is no smooth, everywhere 
spacelike interpolating manifold between the geodesics and Xi Ul 2 ). Of course, the Vaidya-BTZ 
spacetime that we will consider is not static, so the geodesics will not lie on a constant time slice as 
in figure 1. However, even in this more complicated case, we should keep in mind that the “diagonal” 
geodesics (AC) and {BD) will not satisfy the homology condition because the interpolating manifold 
cannot be everywhere spacelike. So, the HMI will be computed as in eq. (1.1) in the rest of the 
paper. 

The motivation for the current work is twofold. Firstly, even though similar models have been 
investigated extensively in the literature (see for example [15]-[22]), the main focus was on non¬ 
compact CFTs, (i.e. dual to Vaidya-BTZ in Poincare coordinates). Significantly less attention has 
been directed to compact CFTs, where finite size effects may well result in unexpected behaviors 
for the HMI, since the generalization to compact space is generically non trivial. For instance, it is 
well known that the HMI in non-compact holographic CFTs typically exhibits a “bump” (see [18], 
[19]), so it is interesting to investigate whether more bumps, or qualitatively different behaviors, 
can appear in the compact case. Secondly, there is a sharp contrast between this, typically non 
monotonic, evolution and the evolution in the adiabatic approximation of the quench. As will also 
be explained later, the adiabatic approximation treats the system as always being in a thermal 
state with the temperature slowly increasing. The dual gravitational description of this process is 
given by a BTZ black hole with its radius slowly increasing (since the radius r+ is proportional 
to the temperature T, T ^ r_|_). We intuitively understand that, in this case, the “growth” of the 
black hole radius will “push back” the boundary anchored spacelike geodesics (since they cannot 
penetrate the horizon of a static black hole, see [23]). Combined with the fact that longer geodesics 
reach deeper into the bulk, this suggests that non monotonic behaviors for the HMI cannot appear 
in the adiabatic approximation. For example, we can see that if the “connected” configuration of 
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geodesics dominates (i.e. / = 0) initially, it will remain so during the entire evolution, in contrast 
to the behaviors observed in [18] and [19]. It is interesting to investigate this problem in order to 
understand why these differences occur from a geometrical point of view. 

For convenience, we are going to set the AdS radius of curvature LAdS to LAdS = !> also 
ignore factors of 4GAr below, loosely identifying the length of spacelike geodesics with Si. So we 
should keep in mind that the HRT formula for the HEE is a good approximation when quantum 
and higher curvature corrections are suppressed, since it is an order 0{G^) result in classical 
Einstein-Hilbert spacetime (with Gm small). In other words, recalling that the central charge c of 
the boundary GET is given by the “Brown-Henneaux relation” c = , we will be working in 

order 0{c), in the limit of large c. As usual, the speed of light c and Planck’s constant h are set to 
1, c = h = 1. 

The plan of this paper is as follows. In section 2 we start by presenting the Vaidya-BTZ 
spacetime and its spacelike geodesic equations. Then, we present the solutions of these equations, 
also describing the geodesic structure of the spacetime in some detail. In section 3 we plot our 
results for the holographic entanglement entropy and the holographic mutual information, before 
turning to the comparison with the adiabatic approximation. In the final section we conclude and 
we discuss some open questions. The appendix contains explicit calculations not presented in the 
main body of the work. 


2 Geodesics in Vaidya-BTZ 


2.1 Vaidya-BTZ spacetime in global coordinates 

As explained in the introduction, we will consider the 2 + 1 dimensional, asymptotically AdS 
spacetime described by the Vaidya-BTZ metric: 

= —/(u, r)dv^ + 2dvdr + r'^d(p^, (2.1) 


where r is the usual radial coordinate with range [0, oo), (j) € [0, 27r) is the angular coordinate with 
period 27r, v € (—oo,+oo) is the “ingoing time” and 


f(v, r) = — m{y) = + 1 — 6{v){r'^ + 1) 


r'^ + 1 = fiir), v<0 

-rl = fo{r), v>0 


( 2 . 2 ) 


where 0{v) is the Heaviside step function. It is clear that /(u, r) interpolates between pure AdSa 
for u < 0 and the BTZ black hole with an event horizon of radius r+ for v > 0, while v = 0 
describes an incoming null shell which collapses and forms a curvature singularity at r = 0 for 
V > 0. More generally, we could consider a shell with non-zero thickness, by replacing 0(v) with a 
smooth function interpolating between 0 and 1, over a time extent characterized by a parameter vt, 
such as ^(tanhv/vt + 1). The Ut 3> 1 regime corresponds to a slow thermalization, which we expect 
to be well approximated by an adiabatic process (see for example [9]). We will focus on the thin 
shell limit (vt —>■ 0), which describes a fast, non-adiabatic thermalization and allows us to obtain 
analytic results, returning to comment on the differences between the two descriptions in the final 
subsection. 

Performing the change of coordinates 

t = v — p+^, p = tan~^r (2.3) 

we can easily see that t is identified with the static time coordinate inside the shell (i.e. for 
u < 0), but not with the static time to outside the shell, which makes sense since the Vaidya-BTZ 
spacetime is not static. In these new coordinates, infalling null geodesics always form angles of 7r/4 
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m. p — t diagrams, whereas this is true for outgoing null geodesics only in the AdS part u < 0. Also, 
t coincides with v on the conformal boundary x K at r —^ oo, where it provides a natural time 
coordinate for the dual field theory and will be denoted by too- Then (2.3) implies that the shell 
starts imploding when too = 0. 

More details about the structure of this spacetime can be found in [24]. 


2.2 Geodesic equations 

Following [24] and exploiting the symmetries of the problem, we can obtain the following first order 
equations of motion for spacelike geodesics: 


L = 

(2.4a) 

E = fa{r)v - r, 

(2.4b) 

r^=E^- (^-l)/a(r), 

(2.4c) 


where a = i,o and the dots denote differentiation with respect to an affine parameter s along the 
geodesics. In the above, L is the angular momentum originating from the spherical symmetry, and 
E is the energy originating from time translation invariance in u < 0 and u > 0. It can be seen that, 
while L is kept constant along each spacelike geodesic, E is constant in u < 0 and u > 0 separately, 
but is discontinuous at the point where the geodesic crosses the shell, i.e. at u = 0. In fact, the 
discontinuity can be computed: 

AA = i?|„=o+ ~ = -{f\y=o+ — /|„=o-)'^|u=o = (2-5) 

with V being continuous everywhere. 

From the second order equation of motion for v 


ldf.2 


( 2 . 6 ) 


we observe that i;|i,=o > 0, which implies that v has one (global) minimum along the geodesic. This 
will be denoted by Vmin- Now, noting from eq. (2.4b) that flipping the sign of E is equivalent 
to the reparametrization s —>■ —s, we will only consider geodesics with E\y^^^ > 0 without loss 
of generality. Similarly, (2.4a) implies that flipping the sign of L is equivalent to reversing the 
direction of the increase of tp, so we will only consider L > 0. Then, eq. (2.4b) and (2.4c) imply 
that r|„ . = L. 

We are interested in computing the entanglement entropy of boundary theory intervals, the 
endpoints of which lie at the same boundary time too- In our case, the HRT prescription dictates 
that this is given by the length of spacelike geodesics with both endpoints anchored on the boundary, 
at the same time too (termed “ETEBA geodesics” in [24]). Generally, there are 3 types of ETEBA 
geodesics in Vaidya-BTZ: 


• Geodesics lying entirely in the AdS part of the spacetime (i.e. in u < 0). These must 
necessarily end at a time slice too < 0 on the boundary. 

• Geodesics lying entirely in the BTZ part of the spacetime (i.e. in u > 0). These can only end 
at too > 0. 

• Geodesics that lie in the BTZ part near the boundary, but reach inside the shell. These, too, 
can only end at too > 0. 
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Geodesics with E\y^.^ = 0 have an additional symmetry under s —> —s, so they necessarily end at 
the same time slice on the boundary (if they do not end at the singularity). In principle we could 
also have asymmetric ETEBA geodesics, i.e. ETEBA geodesics with E\y^.^ > 0, but there are good 
arguments against their existence (see [24]), so from this point on we will assume E\y^.^ = 0. Then, 
geodesics which lie entirely in the AdS or the BTZ part of the spacetime have E = 0 everywhere, 
whereas, adopting the convention that s increases in the direction towards the boundary, (2.5) 
implies that geodesics which cross the shell have E = 0 inside the shell and E < 0 outside. 


2.3 Solutions of the geodesic equations 

Solving eq. (2.4) is straightforward, though somewhat tedious. The details of the calculations, as 
well as the explicit forms of the geodesics, can be found in appendix A; here we just quote the most 
important results. In the following, r^o is the (large) radial cut-off. 

• For too < 0 the geometry is pure AdSa and we obtain the well-known result that, for an 
interval I = [—V^oojV’oo] on the boundary, where ipoo € [0,7r), there is a unique spacelike 
geodesic of length 

= 21n(2roo)-k21n(sinvJoo)- (2.7) 


Before moving on to the next cases, we briefly recall what happens in pure BTZ spacetime. As is also 
explained in [10], fixing the interval X = [—v^coji^oo] on the boundary, the spacelike geodesics can 
wind around the BTZ black hole an arbitrary number of times. This results from the existence of 
solutions that have angular extent greater than 27r, but are anchored at the same endpoints dX on the 
boundary due to the compactness of tp (see figure 1 (right)). More precisely, geodesics with angular 
extent = 2(/?oo + 27rfc, k G as well as those with angular extent = —2ipao + 27r(fc -I- 1), 
all end on dX. However, in the present case the homology constraint is non-trivial; this implies that 
only is homologous to X and is homologous to the complementary interval X. The lengths 


of the geodesics are given by ~ 21n(2roo) -I- 2 In I 

According to the RT prescription, the entanglement entropy of X is given by the area of the bulk 
minimal surface anchored on dX and homologous to X. In our case, candidate bulk minimal surfaces 
include the connected geodesics described above and disconnected surfaces consisting of two pieces: 
the geodesics and the bifurcation surface of the black hole event horizon %. In the latter case, a 
geodesic homologous to X, together with XL, form an extremal surface which is homologous to X. It 
can be seen that there is a value ip^^ > tt 12 oi ipoo separating two regimes: for ipao < Xx HEE 
is given by the length geodesic with angular extent = 2pao, while for p^o > Xx 

it is given by the sum of the length i^g^z g of the geodesic with angular extent = 27r — 2poo 
and the bifurcation surface area areaLXL). This regime was termed “entanglement plateau” and was 
explored in [10].^ The fact that generally Sg Sg signifies that the BTZ spacetime is dual to a 
mixed state in the boundary theory (in particular, to a thermal state). 


sinh(r+:/j*^^^)^. 


Returning to the study of the Vaidya-BTZ case, we firstly note that, due to the fact that the 
spacetime is topologically trivial, the homology constraint trivializes (at least for the analogues of 
the k = 0 geodesics described above^). Thus, the disconnected surfaces are absent in our case and 
we will always have Sg = Sg (signifying that the state is always pure), because dX = dX and the 
geodesics giving the HEEs will be homologous both to X and X. So, from now on the discussion 
will be focused on a single interval X. 


^See also [11] for related work. 

■^More evidence for why A: yf 0 geodesics don’t come into play will be given in subsection 3.1 



Provided that too > <Poo^ (using an obvious notation in analogy with the above), we find 


geodesics which lie entirely in the BTZ part u > 0®. Their lengths are given by: 

^^BTZ = 21n(2roo) + 2 In sinh(r+v5^^)^ . 


( 2 . 8 ) 


• Finally, for too > 0, there exist ETEBA geodesics that cross the shell. Using the radial 
coordinate at which they encounter the shell rg (note that they are symmetric, so there is 
only one such parameter) as a parameter, we find that these geodesics have: 


VJoo(r+,rs,too) = + ipl^{r+,rs,t^), 

where is the angular extent inside the shell 

2(r+ - rgT) 


(2.9a) 


'PL(^+UsUoo) = tan ^ 


- 4r+(r^ - l)rgT + (1 +?’+(2 +'ir‘l))T‘^ ^ 

having set T = tanh(r+too)j and 

I ( J-^r\ - 4r+(r^ - l)rgT + (1 + + r^(2 + 4r2))r2 

= —tanh"^ ' 


(2.9b) 


r+ 


1 - + 2r+rsT 


is the angular extent of the geodesic outside the shell. 

The length is given by the expression 

^(^+ 5 ; too ) (^+ ) j too ) 4“ ^ (^+ 5 1 ^Oo) ^ 

with 

r(r+,rg, too) = In 

giving the length of the AdS part and 


K + 1)T 

(rl + 1)T + 4rs(rsT - r+) 


£°(r+,rg,too) = 21n(2roo) + 21n 


(4 + l)r + 2rg(rgr-r+) 
r+(rl + l)^/(l^ ^ 


(2.9c) 

(2.10a) 

(2.10b) 

(2.10c) 


giving the length of the BTZ part. 

Here we should also mention the following constraints, which we obtain by demanding that 
the geodesics end on the boundary (and not at the singularity): 


- 1 + (r^ + l)\/l - 
2r+T 

In particular, the first inequality in (2.11) implies that 

ri — 1 

r.g > 


< rs < Y' 


2r. 


( 2 . 11 ) 


( 2 . 12 ) 


since 0 < T < 1. The second inequality in (2.11) trivializes when Vg < r+®. 


^Note here that the same condition was found in the Poincare case (see for example [19]), though only for the 
solution with A: = 0, since no geodesics wind around the black hole horizon in that case. 

®These constraints may be thought of as the analogue of the “critical extremal surface” of [21],[22]. They describe 
the allowed parameter space for the geodesics to be anchored on the boundary, and also give a significant amount of 
information about the behavior of certain physical quantities, as will be seen below. 
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From eq. (2.7), (2.8) and (2.10c), we observe that the divergence in the length of the geodesics 
appears as the term 21n(2roo), in which, strictly, we should take the cut-off to infinity r^o —t oo. 
We will invoke the simplest regularization procedure, which involves subtracting this term from the 
above equations. From now on, when we refer to the length of the geodesics, we will always mean 
the regularized length. 

Poincare limit: A good consistency check would be to show that our solutions (2.9) and (2.10) 
reproduce the known solutions in Poincare coordinates, which were first found in [16]. To go from 
global to Poincare coordinates, effectively we need to “zoom in” on the boundary and at the same 
time make the black hole large (so that it does not disappear). That way, the boundary will look 
like and the horizon will look planar. More precisely, we need to scale 

r -)► Ar (=> r+ Ar+, r^o -t Aroo), 

tp^A-V, (2-13) 

V —>■ A“^z; (<;4> t —>■ X~^t) 

and then take the limit A —?► oo. This also has the effect of uncompactifying (p, since the period 
becomes infinite. Then, the Vaidya-BTZ metric (2.1) takes the form 

ds^ = —(r^ — 0(u)r^)c?u^ -I- 2dvdr + r^dx^, (2-14) 


in which we have conventionally denoted the uncompactified direction on the boundary by x. 
In this limit, (2.9) becomes 


‘P^{r+,rs,too) = 


2(r+ - TsT) 


r+-\/Tr+(-4rs -|- r+T -|- dr^T) r+ 


J_tanh-i ( \/^^+(-4rs -f r+T + 4rgT) 


-r+ -h 2rsT 


(2.15) 


and (2.10) becomes 


^^{r+,rs,tao) = In 


/ 4r^r(4T + 2r,(r,r-r+))2 \ 

\TX{l-T^){rlT + ^rs{rsT-r+))) ' 


(2.16) 


We can now see that the solutions obtained in [16] are equivalent to (2.15) and (2.16) after making 
the appropriate redefinitions. 


Graphical representation: In principle, we should now invert the expression for the angle ipaoi 
(2.9), in order to find (r_|_, , tco) for a fixed Pao-, and then replace it in the expression for the 

length •^(r+, Tg, too)) eq. (2.10). This way we can find the evolution of the geodesic length for a fixed 
interval on the boundary, i.e. £(r_|_, i^oo, too)- However, (2.9) cannot be solved analytically for rg, so 
we will solve it numerically in the next section. We end this section by examining some qualitative 
features of the geodesics in global Vaidya-BTZ. 

In figure 2 we have made the contour plot of Pc>o{f+,rs,tao) for a black hole of radius r_|_ = 3. 
The plot continues to rg —>■ oo and too —>■ oo in an obvious way. The coloured curves in the rg — too 
plane are curves of constant poo- We have also emphasized the boundaries of the plot, which can 
be read off from constraints (2.11) for general r+: the thick black line is given by 


too = J^i(r+,rg) 


—tanh ^ 

'r+ 


2r+((r 


+ - + (^+ + l)\/^s + 1) 

l + rX + 2rl{2rl + l) 


(2.17) 


along which p^o = f, and the thin black line is given by 


1 / T* I 

too = '^ 2 (r+,rs) = —tanh“^ 

r+ 


r. 


(2.18) 










5 

- 

8 

- 

24 

______ 2 

2 

_ 11 7T 

24 

_ 3tt 

8 

_ 7T 

5 

— f<x>=vi(r+=3,rs) 

— t<»=fj{r*=3,rs) 

— (<x>=v'2(^=3,rs) 


Figure 2. Density plot of (^oo(’'+ = 3,rs,too)- The continuous coloured curves are curves of constant tpoo, 
as can be seen in the legend. The dotted curves “go around” the singularity, but they have the same spatial 
endpoints as the corresponding continuous ones (on time slices specified by too). The dotted-dashed black 
line, which corresponds to half of the spacetime boundary, can be thought of as a limiting case, with two 
branches merging into one. The boundaries of the plot are the thick black line too = = 3, r^) and the 

thin black line too = i' 2 {r+ = 3,rs). The curve 122 ( 3 , Vs) is a local minimum of ifioa with too constant, while 
= 3, Ts) is a local maximum. 


marking the time after which the geodesics belonging to a family with a particular angular extent 
lie entirely in the BTZ part of the spacetime. 

Let us now focus on a fixed interval of length 2{pca with ipoo < § on the boundary. We observe 
that, for too 0> we get —)■ cx) as expected. Then, as time goes by, gets smaller, reaches a 
minimum, and grows again, until it reaches the boundary i' 2 (r^,rs) of our plot, at which point the 
geodesic just touches the shell. After that time, the geodesics belonging to that particular family lie 
entirely in the BTZ part of the spacetime. In figure 3 (left) we can see this evolution for (p^o = 

It is apparent from figure 2 that, at a specific time, a new family of geodesics will come into play. 
It will have the same endpoints as the family of angular extent 2(/?ooj but it will have an angular 
extent of 27r — 2(pao, he. it will “go around” the black hole. This new family has two branches of 
geodesics, separated by the curve 

too = /r(r+,rs) = ^tanh“^ (^) ’ (2-19) 


namely the right one, with an evolution and BTZ saturation as the one discussed before, and the 
left one, which tends asymptotically to the boundary curve z/i(r+,rs). A graphical representation 
of the evolution of this family can be seen in figure 3 (right), again for (poc = Note that, for 
large enough intervals, the second family appears before the first one has saturated to its BTZ form. 

Finally, it is interesting to observe that the geodesics of the left branch accumulate quickly to the 

-1 

boundary curve J^i(r+,rs) regardless of the interval size and that 121 itself tends to the value 
for large time too- 
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Figure 3. In these plots, the time coordinate t runs upwards, and we have compactified the infinite radial 
direction by plotting p = tan“^ r. The outer cylinder is the boundary of Vaidya-BTZ, the (null) cone is 
the incoming shell and the inner cylinder is the event horizon of the black hole formed after the collapse. 
Here we have chosen r+ = 3 and we are only plotting the part with too > 0. Note once again that t ^ to, 
and that is why the geodesics lying entirely in the BTZ part of the spacetime do not lie on a constant-t 
slice. The green lines are reference lines of ipoo = tt on the boundary (showing that the plots are rotated 
with respect to each other). The red lines mark the endpoints of the boundary intervals in question, i.e. 
(foo = Left: Geodesics belonging to the family with angular extent Right: Geodesics belonging 
to the family with angular extent 27r — The blue geodesic appears first, which then splits into the “left 
branch” (black) and the “right branch” (purple), as shown in the plot. 


For geodesics connecting antipodal points on the boundary, the splitting into radial and non- 
radial ones is just the degenerate case in which the right branch of the k = 1 family coincides 

with a part of the k = 0 family. This happens at t\cusp = —ta-nh” ( —— j , when zzi(r_|_, r^) 

meets /i(r_|_,rs), at which point the curve r’s(tco) for constant = '^/‘2 is not smooth. We thus 
understand that there is a spontaneous breaking of the Z 2 reflection symmetry. 

In a similar fashion as above, higher-/c families, containing geodesics which wind around the 
black hole multiple times, will come into play as too increases. 

A slight difference that occurs for radii r+ < 1 is that the curve i'i{r+,rs) ends on the too 
axis at time ^tanh”^ ^ r^’+i ) (which is always greater than t\cusp)- This means that the family of 
radial geodesics, as well as those belonging to the “left branch” of the higher k families, eventually 
cease to exist (see also the discussion in [24]). 

3 Holographic entanglement entropy and mutual information 

3.1 Holographic entanglement entropy 

The holographic entanglement entropy (HEE) of a boundary interval I at time too is given by the 
length of the minimal ETEBA geodesic anchored on dl at that specific time. In principle, there is 
an infinite number of geodesics we should take into account, since each family (which is labelled by 
k G Z+, see eq. (2.8)), generically contains 2 candidate geodesics, (see figure 3 (right)), apart from 
the family with k = 0 (see figure 3 (left)). We expect that geodesics which wind around the black 
hole many times are longer than those which do not, and, although it is hard to prove this directly, 
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eq. (2.9), (2.10) and (2.8) give strong indications that this is the case. We postpone this discussion 
momentarily, now turning our attention to the /c = 0 families which should give the HEE. 

We can numerically invert (2.9) to find rs(r+, <foo,too) and replace in (2.10) to express the length 
as a function ^(r+, foctao)- The resulting plots of the HEE for various interval sizes are presented 
in figure 4, for r+ = 5 (top left), r+ = 1 (top right) and r+ = 1/4 (bottom). We observe that there 
is a smooth, monotonic interpolation between the AdS value and the BTZ value, always increasing 
since eq. (2.7) and (2.8) imply that £i < to for all ipao- As expected from (2.7), the initial value of 
the HEE for intervals of length 7r/2 is 0, while it goes to —oo as the length goes to 0'. Importantly, 
the linear growth regime initially predicted from general GET arguments in [25], and verified from 
a holographic viewpoint in [14]-[17], [20]-[22] and [26] (among others), is absent in our case. This is 
not surprising: the derivation of [25] relied crucially on defining the 1+1 dimensional quantum field 
theory on as well as starting from a gapped theory before the quantum quench, and ending 
with a GET after the quench (but it was independent of the coupling). Physically, this means that 
initially there were only short range correlations at a scale set by the inverse mass gap, inside which 
the proposed quasiparticles that “carry the entanglement” are produced, and then travel freely to 
infinity with a velocity v = 1. This toy model explains the linear growth and saturation found in 
that case. In the context of AdS/GFT, initially (i.e. for too < 0), the boundary theory state is the 
vacuum of a 1+1 dimensional GET, since it is dual to pure AdSa, and so we start with long range 
correlations. Despite this difference, the linear growth regime can also be derived from the gravity 
perspective by using the separation of scales r+(poo ^ r+too ^ 1 (see for example [22]). The effects 
of this separation of scales will start appearing as r+ becomes larger, and indeed it can be observed 
in figure 4 that a linear regime starts forming when r+ = 5. 

Although the free quasiparticle model is adequate to capture the growth of the entanglement 
entropy of one interval after a global quench, it fails when the subsystem consists of two or more 
disconnected components, as was noted in [27] : it predicts a non monotonic evolution, while the 
holographic prescription and large c GET calculations show that it increases until it reaches its 
equilibrium value. The same contradiction arises for our compact system: the periodic, sawtooth 
evolution predicted by the above toy model is in sharp contrast to the plots in figure 4. We 
thus expect that figure 4 captures the essential characteristics of global quenches in compact 1+1 
dimensional quantum systems, at least for large central charge c®. 

A simple calculation of the behavior of the HEE at early and late times (see appendix B) 
confirms the findings of [29] : for too 0 we find 

+ 1 

£{(poo,r+,too) = 21n(sin(v3oo)) + + C>(iL)> (3-1) 

while for too —t +oo 

^(+oo,r-+,too) = 21n sinh(r+(/Joo)) _ _ too)^^^ 

V»'+ / 3y+/ (3.2) 

+ CI((+oo ~ too)"^)- 

k > 0 families: Returning to figure 2, we make the following observations which hold for general 
r+: Firstly, eq. (2.8) shows that t^ is an increasing function of k. Secondly, it can be checked that 
the qualitative features of +oo(^'++s, ^oo) agree with those of £(r_|_, rg, too); in particular, the curve 
V 2 is a (local) minimum and /r is a maximum (for constant r^), while the parameter space is of 

^Note that we are only considering the constant part of the HEE, having ignored the term 21og(2roo) which 
diverges quicker. 

®Indeed, it was recently shown in [28] that a small enough central charge c leads to an HEE profile in agreement 
with the quasiparticle model of [25], while in the holographic regime (c ^ 1, strong coupling), the profile agrees with 
the one predicted using the “HRT prescription”. 
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0.5 1.0 1.5 , 

^oo 

Figure 4. Holographic entanglement entropy of intervals of total length (from bottom to top) 7r/6, tt/S, 
7r/2, 27r/3, Srr/b, tt, for the cases r+ = 5 (top left), r+ = 1 (top right) and r_|. = 1/4 (bottom), as a function 
of the time on the boundary too- 


course bounded by vi and V 2 - It can be shown by a direct calculation that the length I increases 
along curves rs{too) of constant ipoo, as well as along vi, and /r (see appendix B). In fact, along 
V 2 the length is 

£(r+, too) = 2 In sinh(r+too)^ (3.3) 

and along vi the length is 


€(r+,too) = 21n 


-1 + + (1 + r^) cosh(r+too) 


2r+to 


(3.4) 


for r_|_too ^ 1. It is interesting to note here that in the r+ < 1 case, after j/i(r_|_,rs) terminates, the 
length along the rs = 0 axis is also given by eq. (3.3) above. 

In principle, it is possible that for some too, the geodesics belonging to the k = 0 and k = 1 
families exchange dominance, i.e. the minimal geodesics belong to the k = 1 family. However, the 
allowed parameter space for this to happen can be significantly restricted by the above observations. 
Firstly, note that we can show that the length increases along the k = 0 families of fixed tpoo, as well 
as along the left and the right branch of fc > 1 families (see appendix B). Also, since for constant 
too the length is monotonically decreasing from ^ to 122 , we can see that when they first appear, 
the higher k families will never be shorter than those belonging to lower k families, and this will 
remain so for the right branch geodesics. This confirms the intuition that geodesics winding around 
the black hole are always longer than those which don’t. Regarding the left branch, which contains 
geodesics reaching deep inside the shell, we observe from figure 2 that they accumulate quickly near 
the boundary curve vi, which contains geodesics with lengths growing almost linearly with time too 
(see eq. (3.4)). So we expect them to be even longer, although we cannot rule out the possibility 
that something counterintuitive happens at the beginning of the left branch. Similar arguments 
can be applied to families with general fc > 1®. 

®Note also that in some cases the geodesics may even be excluded from the possibility of giving the HEE due to 
the homology constraint. 
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Figure 5. Time evolution of holographic entanglement entropy of fc = 0 (blue), and the left branch (black) 
and right branch (black, dashed) of fc = 1 families of geodesics. The radius of the black hole is r+ = 5 and 
the intervals on the boundary are of total length: 7r/2 (top left), Stt/G (top right), and tt (bottom). 


After an extensive numerical search of the parameter space, the general behavior is found to 
agree with our intuition, and can be seen in figure 5. Thus, the saturation to equilibrium is found 
to be continuous. 


Rate of growth: An interesting quantity to compute is the dimensionless rate of growth 


^{too) 


1 d^T 
dtoo 


(3.5) 


introduced in [21] and [22] for a boundary region I in a general d-dimensional boundary. In the 
above, A is the entangling surface of I, i.e. A = area{dCL) and Seq is the equilibrium entropy 
density. Seq is defined such that the difference of the equilibrium entropy from the vacuum entropy 
A5'eg takes the form /S.Seq ~ SeqVoliT) for regions T with large effective size Reff compared to the 
equilibrium scale l/r+ set by the radius of the black hole. The definition of the rate of growth 
was motivated from the fact that, in a very large class of Vaidya-type geometries with non-compact 
boundary, there exists a universal “post-local-equilibration linear growth” of the form 


Sx{t 

oo ) = VESeqAt 

C)0 ■ ■ ■ 5 (3.6) 

where is a shape-independent dimensionless quantity. This linear growth regime appears when 
Reff 3> too When too 1/?'+ a quadratic “pre-local-equilibration growth” was observed. 

Specializing to the Vaidya-BTZ case in Poincare coordinates and translating to our notation, 
we get: Seq = A = 2 (reflecting the fact that a boundary interval has 2 endpoints) and ve = 1 
(a universal result for the linear growth regime in non-compact 1+1 dimensional CFTs). So we get 

^ 1 (3-’^a) 

At = 2r+Y, r+:poo > T’+too > 1 (3.7b) 

our analysis we consistently ignore the factor of 1/AGjsi-, as explained in the introduction. 
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and 




1 d£ 


2r4 


(3.8) 


+ d^oo 

We also note that extensivity for the equilibrium entropy is derived by the expression for the length 
of a spacelike geodesic in BTZ spacetime: eq. (2.8) implies tsTZ ~ ^+(2:poo) when r^Lpoo ^ 1- 
Natural arguments led the authors to suggest that 93 < 1 for the class of systems they studied 
and in appendix C it is shown that this holds in the Vaidya-BTZ case in Poincare coordinates. 
However, it is easy to see that we can have 93*-^^ > 1 in the case of global coordinates for small 
enough For example, we can compute the rate of growth of a boundary interval with ip^o = 7’'/2 
at the time when the geodesics stop being radial: the time derivative of eq. (3.4) evaluated at time 
I 1 ^ f '^+ \/^"^+2 \ 

t\cusp = T^tanh ( ) gives 


93(2)|,„,p = 


2{rl + l)-^r 




+ 2 


(3.9) 


This implies that 93^^^|cMsp > 1 for r_|_ < \/— 1 < 1, while 93*'^^|cusp 1 from below quickly as 

r+ grows larger. 

This does not imply any violation of causality, but only that the definition (3.8) ceases to have 
a nice interpretation as the rate of growth of HEE when r+ is small enough. More precisely, a small 
r+ cannot be interpreted as the entropy density Seq if our system is compact, since extensivity will 
fail even for the biggest interval that we may consider (we necessarily have tpoo ^ Indeed, we 

see that the divergence at r_|. = 0 comes from the normalization l/2r+. This also makes sense if we 
recall that the extensive behavior of the equilibrium HEE comes from the part of the geodesic which 
“lies along the horizon”, so a small radius black hole cannot have this effect in a compact boundary 
spacetime. On the other hand, in the flat boundary case, r_|_ can always be interpreted as the 
entropy density s^q, since we can consider arbitrarily large boundary intervals, able to compensate 
for the small . Geometrically, a spacelike geodesic anchored on a large enough boundary interval 
will always contain a part lying along a flat, non-compact black hole horizon. 

3.2 Holographic mutual information 

We now turn to the evolution of the holographic mutual information, where we can observe im¬ 
portant deviations from the adiabatic approximation and a generalization of the previous findings 
in the Poincare case ([18] and [19]). All qualitatively different behaviors are presented in figure 6 
below, for r+ = 5 (without loss of generality). 

We observe that the evolution of the HMI strongly resembles the flat boundary case in the 
upper line and the left panel of the middle line in figure 6. More specifically, considering only same 
length intervals and varying their separation, the classification of [19] includes 4 different behaviors: 

• the HMI is 0 at all times 

• the HMI starts from 0, becomes positive for some time and ends at 0 

• the HMI starts from a positive value, forms a “bump” and ends at 0 

• the HMI starts from a positive value, forms a “bump” and ends at another positive value, 
smaller than the initial. 


^^Numerical calculation of for various black hole radii r+ and boundary interval sizes tpoo is in agreement 

with the following arguments. 
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Figure 6. Time evolution of holographic mutual information for = 5. Upper line: HMI of two 
intervals of lengths 7r/2 and 7r/2 (left), and Stt/S and Stt/S (right), separated by intervals of length 7r/24 
(black), 7r/8 (purple), 77r/24 (red), 57r/12 (yellow) and 7r/2 (blue). Middle line (left): HMI of two 
intervals of lengths 77r/16, separated by intervals of length 7r/24 (black), 7r/6 (purple), 77r/24 (red), 7r/2 
(yellow) and 77r/12 (blue). Middle line (right) and lower line: HMI of two intervals of lengths 97r/16 
and 97r/16 (middle line right), 77r/16 and IItt/IG (lower line left), and Stt/IG and IStt/IG (lower line right), 
separated by intervals of length 7r/24 (black), vr/G (purple), 77r/24 (red), Stt/S (yellow) and 77r/lG (blue). 
Note that some of the above lines are always 0. 


Making the two intervals unequal has the effect of “flattening” the bump, but no other interesting 
phenomena arise. 

These 4 cases appear in our study of the compact boundary as well, when the sum of the two 
intervals is less than half of the total boundary. For instance, the middle left panel depicts the HMI 
of two equal intervals of length 2(/)oo = Ttt/IG, the separation of which increases from top to bottom 
(for a large enough separation the HMI is always vanishing). This holds until both intervals are of 
length 7r/2 (see the upper left panel). Making the two intervals unequal but keeping the same total 
length produces a deformation similar to the Poincare case (see the upper right panel). 

Finite size effects seem to take over when the sum of the lengths of the two boundary intervals 

[30] it was shown that the HMI of strip regions in higher dimensions obeys the same classification. 
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exceeds tt. When this sum is large enough, the final value of the HMI is greater than the initial, as 
can be seen in the middle right panel and the lower line in figure 6 (which show the evolution of 
the HMI for intervals of total length Qtt/S, with their separation increasing from top to bottom in 
each panel). Of course, the initial and final values are dictated by the equilibrium formulae (2.7) 
and (2.8), so it is more interesting to observe the interpolation between them. 

For two intervals of length 97r/16, we observe that the bump is present when their separation 
is small enough, but at maximum separation we see a smooth monotonic increase. This is a result 
of the effective degeneracy of the parameter space due to the symmetry of our configuration, and 
it holds until one of the complementary intervals has the same length as one of the intervals in 
question (see the red line in the middle right panel). After this point the bump starts to appear 
again. 

The lower line in figure 6 shows how the above considerations are altered when the two inter¬ 
vals have different lengths. In particular, from the lower right panel we can understand that the 
monotonic behavior appears when one of the complementary intervals becomes equal to one of the 
intervals in question. When their separation is small, we see a quite deformed bump. When their 
separation is large, the HMI decreases in the beginning and then increases up to the saturation 
value. When their separation is large enough, the HMI can even reach 0 and remain there for some 
time before it starts increasing. For even more asymmetric intervals, the HMI can begin from 0 
and then start increasing at a specific time^^. 

Many of the above observations result from the compactness of the boundary. In order to gain 
a clearer understanding of how this happens, some of the curves of figure 6 (plus one not included) 
are presented in figure 7, together with a graphical representation of the boundary intervals in 
question. 

The behaviors presented in the top line in figure 7 can be understood as “complementary” 
to some of the behaviors shown in the middle left panel in figure 6 (the purple and the red line, 
qualitatively). Indeed, having partitioned our system into four intervals, the HMI of two of them 
equals minus the HMI of the complementary intervals (always bounded by 0 from below), see eq. 
(1.1). More specifically, the intervals in the middle left panel in figure 6 are of total length 77r/16 
each, whereas in the top line in figure 7 the complementary intervals are of total length 77r/16 each. 
These behaviors are absent in the non-compact case, unless we consider infinite size intervals. 

The behavior shown in the bottom left panel in figure 7 is also a result of the compactness of 
the boundary. Although a bump is present, the biggest complementary interval is slightly bigger 
than each of the two intervals in question, and so the HMI decreases a little before saturating at a 
value larger than its initial one (the fact that only the biggest complementary interval contributes 
to the decrease can be seen from the vertical dotted lines, which mark the saturation times of the 
various intervals). In contrast, in the non-compact boundary case, the corresponding interval would 
be larger than the sum of the intervals in question, resulting in a decrease of the HMI to a value 
always smaller than its initial one. 

Now, the evolution shown in the bottom right panel in figure 7 is interesting because no bump 
or “inverse bump” is present. As explained above, it can be understood as an intermediate be¬ 
havior between the cases in which the contribution of one of the intervals in question cancels the 
contribution of one of the complementary intervals. Two stages of increase are distinguishable in 
this plot, occurring between the dotted lines for the following reason: the size of the small interval 
57r/16 is bigger but comparable to the size of the separation tt/G and the size of the big interval 
IStt/IG is slightly bigger than the size of the separation 177r/24. In the beginning of the evolution 
we can ignore the two biggest intervals because the corresponding geodesics lie mostly in the AdS 
part of the spacetime and get deformed slowly; thus their contributions almost cancel each other. 

^^See bottom right panel of figure 7. 
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Figure 7. Time evolution of holographic mutual information for r+ = 5. Upper line: HMI of two 
intervals of lengths 57r/16 and IStt/IG (left), and Stt/IG and IStt/IG (right), for complementary intervals 
of lengths Ttt/IG in both cases. Lower line: HMI of two intervals of lengths Qtt/IG and Qtt/IG (left), and 
Stt/IG and IStt/IG (right) for complementary intervals of length tt/G and 177r/24 in both cases. In each 
panel, the intervals in question are graphically represented by the thick black parts of the total boundary 
circle. The vertical dashed lines mark the saturation times of the different intervals that come into play in 
computing the HMI (note that in the top line and the bottom left panel two of the dashed lines coincide 
because, in each case, two of the intervals are of equal length). 


So, after the smallest of the two complementary intervals has saturated (first dotted line), the main 
contribution to the HMI comes from the interval of length Stt/IG, leading to the first clear stage of 
increase. In an analogous way we can understand the second stage of increase as the competition 
between the remaining two intervals. We thus see that this is also a direct result of the compactness 
of the boundary. Similar arguments can be applied to most of the curves presented in figures 6 and 
7. 

Finally, from figure 6 we note that for any configuration, as we bring the two intervals close 
together, the HMI seems to obtain a fixed shape, which is then only shifted to larger values. This 
happens because, to first order, the HEE of a small interval is just a (negative) constant (see 
equations (2.7), (2.8) and also figure 4). 

3.3 Comparison with adiabatic approximation 

Finally, let us compare these results with what we would get using the adiabatic approximation. 
From the boundary theory point of view, we slowly inject energy to the system, so that it (ap¬ 
proximately) stays in equilibrium during the whole evolution. Holographically, this amounts to 
effectively making the infalling shell thick enough, so that spacetime looks like a BTZ black hole 
with its radius slowly increasing in time. In other words, when the mass function m(v) in eq. 
(2.2) is slowly increasing^^, we can approximate our fixed-u spacetime slice with a BTZ spacetime 
slice, where = ^m(u). The spacelike geodesics giving the HEE will then have their equilibrium 

^■^Note that the null energy condition does not allow m(v) to decrease. 
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Figure 8. Time evolution of holographic mutual information among two intervals of lengths 137r/24 
separated by an interval of length 27r/24, for r+ = 1. The continuous curve represents the case of an 
infinitely thin shell, while the dashed curve represents the case of a shell with effective thickness Vt = 1- 


forms, i.e. eq. (2.8), with the replacement r_|_ = ^Jm{v). This is equivalent to forcing our system 
to remain in equilibrium. In the rest of this subsection, will refer to the final black hole radius 
Ty = lim„_>oo \/^fn(v) and we will use m(v) to describe the intermediate stages of the evolution. 

Pure AdSa is obtained for m(v) = —1. As m(v} increases, a naked conical singularity forms 
until m(v) = 0, while for m(v) > 0 an event horizon appears^®. It is natural to consider the mass 
profile 

m(v) = — tanh f — j + —, (3.10) 

smoothly interpolating between to = — 1 and to = with vt corresponding to the effective 
“thickness” of the shell. The limit Ut —>■ 0 reproduces the previous analysis of the instantaneous 
quench. 

In [9] the adiabatic approximation was checked against numerical calculations for ry = 1 and 
Vt = 1, and a very good agreement was observed for the HEE of intervals of various sizes. This 
implies that we can compute the HMI using the adiabatic approximation with the mass profile 
m{v) = tanh(r;) and we are guaranteed to obtain trustworthy results for a slow enough quench. So 
we are going to consider this case, with our analysis also persisting to other values of ry. 

Figure 8 confirms the expectations outlined in the introduction, in the specific case of two 
intervals of lengths 137r/24 separated by an interval of length 207r/24: even though the infinitely 
thin shell produces a “bump” in the HMI, the thick shell predicts a smooth monotonic evolution. A 
similar monotonic evolution between the initial and final value of the HMI was also confirmed for 
many other choices of intervals, in contrast to the variety of different behaviors presented in figure 
6. 

In order to highlight the differences between the adiabatic and the non-adiabatic evolution, 
both cases are presented in figure 9. In particular, the top left panel shows all of the geodesics that 
come into play in the non-adiabatic evolution, from too = —7r/4 up to tco = 2 (for clarity, the shell is 
not depicted). The top right panel shows the analogous process using the adiabatic approximation. 
In the bottom panel we can see the projection of the four intermediate geodesics on the p — plane, 
with the color coding indicating the rate of change of the affine parameter s along the geodesic. We 
have imposed a cut-off for the geodesics at r = 10. 

In figure 9 we can observe the geometric differences between the two processes, the sharpest 
being the way the geodesics are “deformed”: in the non-adiabatic evolution we see that the “pertur¬ 
bations” propagate from the endpoints towards the middle part, whereas in the adiabatic evolution 

^®The Hawking-Page transition for 0 < m{v) < 1 does not concern us here. 
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Figure 9. Upper line: The conventions are as in figure 3. Here, too runs from —7r/4 up to 2 and the 
black hole radius is r+ = 1. Neither the shell nor the black hole horizon is depicted. The boundary is 
partitioned in 4 intervals of lengths 207r/24, 137r/24, 137r/24 and 47r/24, as indicated by the dashed lines. 
The corresponding spacelike geodesics are presented for too = —rr/d, 0, 0.5, 1, 1.3 and 2, for the non- 
adiabatic (left) and the adiabatic case (right). Note that the black geodesics seem to be discontinuous only 
due to numerical inaccuracies. Lower line: Combination of p — ifi geodesic projections in the non-adiabatic 
(left part) and the adiabatic case (right part) separated by the dashed line, for an interval of length 207r/24 
at times too = 0, 0.5, 1 and 1.3. The color coding indicates the rate at which the affine parameter s changes 
along the geodesics: the quicker the color changes, the faster s increases. More specifically, s increases by 
2 as we go from red back to red. Note that the geodesics get deformed in an outwards direction as too 
increases, and that the first two geodesics in the non-adiabatic case almost coincide. In all of the above 
plots, we have imposed a cut-off for the geodesics at r = 10. 


the singularity “pushes” the geodesics mainly at their middle part. This produces a tension in the 
following sense: big intervals feel the effect of the singularity formation early in their adiabatic 
evolution, since they reach deep inside the bulk geometry. In contrast, in the non-adiabatic case, 
the length increase is slow in the beginning because the deformation starts from the UV (i.e. the 
bulk asymptotic AdS region). The length increases significantly only when the shell has “collapsed 
enough” so that the geodesic can feel the effect of the mass. This can be seen in figure 9 (lower line), 
taking into account the fact that the next to last geodesics in both cases have the same length. In 
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the non-adiabatic case we observe that the gravitational pull of the shell makes the geodesic long 
in the region close to the shell, compensating for the middle part which is in pure AdS. Similar 
arguments apply to small intervals, too: in the instantaneous quench case the shell passes through 
the geodesics in the beginning of the collapse, and so they quickly saturate to their equilibrium 
form. In the adiabatic case, it takes them a long time to “feel” the effects of the singularity, since 
they are close to the boundary^®. 

It is important to note that making the parameter vt smaller while still employing the adiabatic 
approximation does not result in the HMI interpolating between the behaviors presented in figure 8. 
This shows that the instantaneous quench cannot be modelled by an adiabatic process, as expected 
from the above discussions. However, the fact that in both cases there is only one family of geodesics 
which gives the HHE of a boundary interval suggests that this should also hold in the numerical 
solutions of the cases with 0 < Ut < 1. We expect that, as vt grows from 0 to 1, the cusps at the 
places where a geodesic crosses the infinitely thin shell get “smoothed out” gradually, until they 
look similar to the geodesics in the adiabatic approximation. Accordingly, we expect the HMI to 
interpolate between the dashed and the continuous curves in figure 8, with the bump starting to 
form at some intermediate stage. 

4 Conclusion and discussion 

To sum up, in this paper we have explored the evolution of the holographic entanglement entropy 
and the mutual information after a quantum quench with the dual description of Vaidya-BTZ 
spacetime. Figure 6 presents a classification of the qualitatively different behaviors of the HMI, 
confirming the behaviors that can be intuitively understood by the arguments in subsection 3.2. 
We also highlighted the geometric differences between the instantaneous and the slow quench. 

We finish our work with a few comments and open questions. Firstly, it is interesting to note 
that the recent proposal of [27] for calculating the (vacuum subtracted) entanglement entropy using 
the “entanglement tsunami” picture of [21],[22] seems to be able to capture many qualitative features 
of the HMI even in our compact boundary case. Although some modifications are required in order 
to achieve good agreement (it only describes the linear part of the HEE in the non-compact case, 
which is absent here), it provides strong support for the usefulness of an “entanglement tsunami” 
toy model. The importance of such a description was also hinted in our comparison of the thin 
shell and the thick shell spacetimes: we saw that the “deformation” of the geodesics propagated 
from the endpoints and the midpoint respectively, resulting in the different evolutions of the HMI. 
Now, we expect that a slightly thick shell will give rise to a tsunami with a slightly thick wavefront. 
However, this picture should break down as the thickness of the wavefronts becomes comparable to 
the size of the intervals in question (which are even comparable to the whole system in the compact 
case), and indeed the behavior observed in the adiabatic approximation was quite different from 
what a tsunami picture would suggest. 

Finally, we should mention that there are many generalizations of the simple Vaidya-BTZ model 
that one could consider. Adding charge and angular momentum to the black hole would be the most 
natural generalization, with the geodesic structure of the spacetime becoming complicated enough 
to allow for novel phenomena in the time evolution of the HMI. More importantly, a systematic 
study of higher dimensional Vaidya-type spacetimes in global coordinates, along the lines of [22], 
would be highly desirable. This would illuminate the possibility of observing genuinely finite size 

^®We should mention that the recent work [28] gave strong evidence that the entanglement scrambles maximally 
in the holographic systems considered by the authors, and so the bump in the HMI is only a result of the different 
saturation times of different intervals. Since the bump disappears in the adiabatic approximation in the non-compact 
case as well, it can be inferred that this is what produces the tension between the fast and the slow quench, as was 
also explained from a geometrical point of view above. 
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effects, invisible to the latter analysis. Of course, a direct calculation of these processes with CFT 
techniques, analogous to those recently developed for the non-compact case (see for example [31], 
[32]), would allow us to obtain a spherical perspective on thermalization. 
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A Solving the geodesic equations 


Here we present the solution of the geodesic equations (2.4). First they must be solved inside the 
shell (i.e. for v < 0), then outside (i.e. for v > 0), and finally the solutions must be patched together 
on the shell u = 0, taking care to keep only those that end on the boundary. Most of the following 
results were also obtained in [29] . 


Inside the shell: Recalling that we only consider symmetric geodesics, the energy should vanish 
inside the shell, E — 0, and we can set the affine parameter s to be equal to 0 at the minimum v 
point Vmin- Then, for u < 0, the geodesic equations (2.4) become 


L = 


(D)2 + r 

(’‘‘1^ = - (^ -1) ((C)" + 1). 

The solution of (A.lc) subject to r(s = 0) = L is 


(A. la) 
(A. lb) 

(A.lc) 


r<(.,L) = y (r" + l)™M2») + L»-l 
and the solution of (A.lb) subject to the same condition is 


v\s,L,Vmin) = tan 


-1 


f Eis,L)-L \ 
\1 -I- r’'{s,L)L J 


+ V 


min • 


(A.2) 


(A.3) 


Without loss of generality we can set (p{s = 0) = 0, and then (A. la) implies that 


(/3*(s, L) = tan ^ 


tanh(s) \ 

■ 


(A.4) 


The parameter value sq at which the geodesic encounters the shell u = 0, will actually give us (half 
of) the length of the AdS part v < 0. Setting E{sq^L) = r^, from (A.3) evaluated at sq we get: 


'^min 


= tan ^ 


L-Ts 
1 -|- r^L 


(A.5) 


We can now invert (A.2) and (A.4) to find (half of) the proper length inside the shell E = Sq and 
(flsheii as functions of the two parameters rg and L: 



+ l + 2y/{r’^+ l){r‘^ 
L2 + 1 



(A.6) 
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and 


(/^Ueii = tan ^ — 


1 r?-L2 


(A.7) 


Outside the shell: In region i; > 0, the geodesic equations take the form 

L = 

E = {{r°f - rl)v° - f°, 


{f°Y =E^- 




(A.8a) 
(A.8b) 

(A.8c) 


(r^ +1) 

The relation E = -^—^'y|i;=o (2-5), the continuity of v across the shell and the explicit expression 

for v"^, eq. (A.3), give the energy in the region ?; > 0 

1 ) Ir^s-L^ 


E = - 


{rl 


2rs y ^2 + 1 
Also, evaluating eq. (A.8b) at r; = 0+, we find 

E{2r‘l — rl + 1) 


I shell — 


(4 + 1) 


= /3. 


(A.9) 


(A.IO) 


We observe that the velocity of the geodesic just after it crosses the shell is positive when > 
(r^ — l)/2 and negative when < (r^ — l)/2. Setting now s = 0 at the point where the geodesic 
crosses the shell, i.e. r°(s = 0) = rg and v°{s = 0) = 0, we can solve (A.8c) with the aid of eq. 
(A.9), to get 


1 


r°{s,a, P,r+,rs) = -\/—2a + 2{2r‘l + a) cosh(2s) + d/Ir* sinh(2s). 


(A.ll) 


where we set a = E'^ — for brevity. 

Eq. (A.8c) determines the behavior of {f°Y, and we can thus find the necessary and sufficient 
conditions which the geodesics should obey in order to end on the boundary: 

• If there is no turning point, i.e. (r°)^ > 0 always, the geodesics should emerge from the shell 
with a positive velocity, i.e. r°{s = 0) > 0. Eq. (A.8c) implies that the minimum of (f°Y is 
at y^Lr+, taking the value E^ — {L — r+Y there. This means that when E"^ > [L — r+)^, we 
necessarily have r1 > {r\ — l)/2. 

• If there is a turning point, i.e. {f°Y = 0 at some point, this should be at a value r of r with 
f < Tg, so that the geodesics emerging from the shell with positive velocity don’t get affected, 
while the ones emerging with negative velocity do change their direction. In this case eq. 
(A.8c) implies that f = -y/—+ ''’s- 

Keeping these conditions in mind, we can proceed to solve eq. (A.8a): 

Lr| tanh(s) 


ip°{s,^,L,r+,rs,<Ps) = — tanh 

r+ 


-1 


SsYs + P tanh(s)) 
Similarly, the solution of eq. (A.8b) subject to v°{s = 0) = 0 is 


■ 


shell ’ 


(A. 12) 


u°(s,/3,£',r+,rg) = 


1 

r+ 


tanh 


tanh(s) 


-r\ + rs{rs + ,5tanh(s)) 


- tanh 


-1 ( r+(rg -r°(s)) 
4 - rsr°{s) 


(A.13) 
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We can now take the limit s —> cx) in (A. 13) and solve for the energy E in terms of r+, rg and 
the time on the boundary too = v°{s —>■ oo). Explicitly, 

(rl + l)(r+ - rgT) 

T(2r2 + r^ + l)-2r+rg’ ^ ^ ^ 

where T = tanh(r_|_too). So, through eq. (A.9) we can express L as a function of r_|_, and tco, 
as well as the quantities sq </?|shell- (Half of) the extent of the boundary interval will be given 

by the s —>■ oo limit of (A. 12), while (half of) the proper length in the BTZ part can be found by 
inverting eq. (A.11). Specifically, assuming r^o to be a large radial cut-off, we find 

4r*^ 

"P(^»") = 2,-; + a + 2». 

Expressing all of the above as functions of r+, and too, we are led to the results (2.9) and 
(2.10) presented in section 2.3. In an analogous way, the conditions for the geodesics to end on the 
boundary boil down to eq. (2.11). 


B Early and late time evolution of holographic entanglement entropy 

Here we present the calculations of subsection 3.1. 

Firstly, it is convenient to make the change of variables rg e-)- a defined by: 


- 1 -I- (r^ -I- l)-v/l - tanh(r+foo)^ 


(1 — a) -I- a~ 


r+ 


(B.l) 


2r+tanh(r_|_too) tanh(r+too) 

This way a interpolates between vi for a = 0 and 1^2 for a = 1, while it simplifies the constraints 
(2.11) to 

0<a<l. (B.2) 

So, we can express ipoo and £ as functions of a in the following way: 

a/2(1 - a) sinh(r+too/2) \ 


(Pooir+,a,too) = tan 


1 / 0(2 — a -I- a cosh(r_|_too)) 


1 , -1 / \/2a{2 - a + acosh{r+too)) sinh(r+too/2) 

—tanh - --— --- 

r+ \ 1 — o-I-acosh(r+too) 


(B.3) 


and 


^(r+,a,too) = In 


/ ((—3(a — 1)^ + (—3a^ + 4a + l)r^) sinh | 

(-b) 

1 + ((o - 1)^ + (o^ + l)r^) sinh | 

(2A=) 

l)^\ 

1 8r:^(-l + (2 - a)a(r^ + 1) + (1 + ( 

i(—2 + o + or^)) cosh(r+too)) 




(B.4) 

Fixing r_|_, we are looking for paths a(too) such that ipoo in eq. (B.3) is constant. So, forcing 
the time derivative of eq. (B.3) to vanish, we get the relation 

ar+ coth ^ ^* 2 °° ) (^+ ~ 1 + cosh(r+too) + a(2 -f — 2 cosh(r+too)) + 2a^(r^ -I- 1) sinh ^ ^'*' 2 °° ^ ) 


a = 


{a - l)(r+ -I- 1 - cosh(r+too) + 2a(r^ -I- 1) sinh ) 


(B.5) 

where we denote by prime the derivative with respect to too- Using (B.5) in the time derivative of 
eq. (B.4) we find that along paths of constant (poo 


2r+(r^ -I- 1)(1 - a) sinh(r+too) 


- 1 -I- (r^ -I- 1)(2 - a)a -f cosh(r+too) + (r+ -I- a(-2 -f a(r^ + 1))) cosh(r+too) 


. (B.6) 
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Now, recalling the restriction (B.2), it can be shown that l'\^p^ > 0 for any t^o > 0, with 
the equality holding only when t = 0 or a = 1, i.e. when the geodesics start evolving and when 
they reach their final BTZ form. This shows in particular that the time derivative of the HEE is 
continuous. 

We can see in more detail the early and late time behavior of the HEE by expanding pertur- 
batively a(too) around the points too = 0 and too = <Poo- Specifically, we see that if 


^(^oo) - 


1 


2 ^ 6 + (12 + r^) tant/)^ 


4tan^ 


t + 

''oo ' 


48 tan (j)‘^ 


4 + o{tl), 


then 


<foo{too) = + 

and substituting eq. (B.7) in the length (B.4) and expanding in tea, we find 

£ = 2\n{sin(l)) + +0(1^), 

as was also stated in subsection 3.1. 

Similarly, expanding a (too) in ((i — too we find 


^(too) - 1 


sinh(r+(())^2r+ tanh(r+(()) r ———— 
- - 1 - 


cosh(r_|_0) — 

2sinh(r+<())(-3r2 +(r^ + l)tanh(r+<()))^_^ ^ ^ ^ ^ ^ 3 / 2 , 

3r+cosh(r+</>)-l tool + ^U^ tool b 


so that 


7^oo(too) = ^ + - teaf). 

Substituting in (B.4) and expanding in ((> — too we find 


£ = 21n ( — sinh(r+(())J — 


2y/2{r^ + l)-^tanh(r+(()) 


(B.7) 

(B.8) 

(B.9) 


(B.IO) 


(B.ll) 


(0-too)"/' + O((<(>-too)'). (B.12) 


An analogous analysis can be performed in the Poincare case, resulting in similar expressions 
with some slight changes in the coefficients of the above formulas. They can be computed exactly 
by taking the Poincare limit (2.13)), reproducing the Vaidya-BTZ case of the general results of [22] 
in the early growth and saturation regimes. 


C Rate of growth in Poincare Vaidya-BTZ 

Taking the Poincare limit of the corresponding expressions in B, we obtain the time derivative of 
the length along paths of constant ipea 

2r+(l - a) sinh(r+too) 


e'P\ 


where now 


rs = r+ 


(2 — a)a + 1 + (1 + a^) cosh(r+too)' 


-(1 — a) + a- 


(C.l) 


(C.2) 


2 tanh(r_|_too) tanh(r+too) 

and so a is still bounded by 0 and 1. Now, since the derivative of eq. (C.l) with respect to a at 
tea we have 

__ 1 _ 1 _ / r . i- \ 

< 1. (C.3) 


constant too is always negative, we find that is a decreasing function of a. Thus, at every 


< T^£'^\a=o = tanh 


f r+to 


V 2 


We conclude that, in the fiat boundary case, the rate of growth as defined above is indeed 

bounded by the speed of light. 
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